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We analyze the time profile AT{t) of the temperature difference, measured across a very compressible 
supercritical ^He fluid layer in its convective state. The experiments were done along the critical 
isochore in a Rayleigh-Benard cell after starting the vertical constant heat flow q. For q sufficiently 
well above that needed for the convection onset, the transient AT{t) for a given e = {T — Tc)/Tc, 
with Tc = 3.318K, shows a damped oscillatory profile with period tosc modulating a smooth base 
profile. The smooth profile forms the exponential tail of the transient which tends to the steady-state 
AT(cxd) with a time constant Ttau- The scaled times tosc/to and Ttaa/tn from all the data could be 
collapsed onto two curves as a function of the Rayleigh number over ~ 3.5 decades. Here to is the 
characteristic thermal diffusion time. Furthermore comparisons are made between measurements of 
a third characteristic time tm between the first peak and the first minimum in the AT{t) profile and 
its estimation by Onuki et al. Also comparisons are made between the observed oscillations and the 
2D simulations by Onuki et al. and by Amiroudine and Zappoli. For e < 9 x 10~^ the experiments 
show a crossover to a different transient regime. This new regime, which we briefly describe, is not 
understood at present. 

PACS numbers: 44.25.+f, 464.70.Fx, 47.47.Te 



I. INTRODUCTION 

Many common fluids can be treated as incompress- 
ible, i.e. one can assume that their density, while 
perhaps varying slightly with temperature, is pressure- 
independent. This is one of the assumptions of the so 
called Boussinesq approximation[l]. Throughout this pa- 
per, we refer to a fluid as "compressible" if the non- 
zero value of the isothermal expansion coefficient (3t = 
1/ p{dp/dP)T significantly influences the fluid dynamics 
and thus cannot be neglected. Here p is the density and 
P the pressure. 

The increased interest in understanding the properties 
of a compressible flow is driven by both experimental and 
theoretical work on highly developed turbulence (See the 
review article by L.P. Kadanoff[2] and references therein). 
Compressible supercritical fluids allow to explore flow 
regimes which are not easily accessible experimentally: 
for example, Prandtl numbers Pr of the order of 500 can 
be obtained, while for an ordinary gas Pr is of order 1[1]. 

A standard experimental approach to study flow is to 
use a Rayleigh-Benard (RB) geometry, i.e. a laterally 
confined horizontal layer of fluid filling the space between 
two horizontal surfaces [3-5]. In a typical experiment, 
one applies externally heat with flux density q, and mea- 
sures the temperature difference AT{t) across the layer. 
The observable quantity of interest to theory is the steady 
state value AT(oo) as a function of q. If the fluid prop- 
erties are known, the relations Nu{Ra,Pr) can be ob- 
tained and compared with predictions. Here Nu and Ra 
are the Nusselt and Rayleigh numbers, which represent 
the normalized heat current and temperature difference, 
respectively. 

Generally, a compressible fluid does not obey the 



Boussinesq approximation. However, in a Rayleigh- 
Benard geometry, the Boussinesq approximation is al- 
ways valid provided that the layer height is sufficiently 
small given the value of Pt- And vice versa, even a 
weakly compressible fluid will violate the Boussinesq ap- 
proximation if the layer thickness is large enough [6]. 
Therefore, while the effects due to finite Pt are inter- 
esting on their own, they can also pose a challenge in 
understanding results of experiments focusing on other 
phenomena: for example, studies of highly turbulent flow 
are often performed in tall RB cells, which is desirable 
because of the rapid (~ h^) increase of the Rayleigh num- 
ber with the cell height. 

We are interested in the effects due to non-vanishing 
f3T, and we use supercritical '^He as a test fluid. Super- 
critical fluids provide a unique opportunity to study the 
compressibility-related effects, because /3t exibits a rapid 
divergence at the critical point. As a result, one can con- 
tinuously change the system from a weakly-compressible 
state to a regime where the dynamics is dominated by 
the high compressibility. The analysis presented here ex- 
tends that presented in an earlier paper [7] and focuses 
on the time-dependent transient AT{t). As explained 
there in detail, the apparatus design was chosen so that 
the high compressibility is the only signiflcant deviation 
from the set of assumptions required by the Boussinesq 
approximation. 

Even in the absence of convection, e.g. under zero 
gravity conditions, the dynamics of a highly compress- 
ible fluid can be qualitatively different from an incom- 
pressible system. As was first predicted theoretically by 
Onuki and Ferrell [8] , and then later demonstrated exper- 
imentally by several groups[9], the density-temperature 
dynamics following a small disturbance of the tempera- 
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ture of the walls depends on whether the pressure P or 
the total volume V ( in other words, the space-average 
density < p > ) is kept constant. In the latter case, 
the mechanism called "Piston Effect" dominates the dy- 
namics: Following a small positive step in the wall tem- 
perature, sharp gradients of p and T form quickly at the 
walls and produce "boundary layers" , which expand away 
from the walls into the bulk of the fluid. The "piston" 
time scale for the formation of the boundary layer is typ- 
ically seconds, while the relaxation back to the uniform 
density distribution happen on the diffusive time scale 
which in practice is minutes to hours. These boundary 
layers compress the bulk fluid adiabatically, which gives 
rise to an apparent fast relaxation of the bulk tempera- 
ture: it becomes correlated with the wall temperature on 
the piston time scale instead of that of the diffusion. At 
constant P, or if the comprc;ssibility is small, the density 
field does not respond to the wall temperature variations, 
no boundary layers are formed and the temperature un- 
dergoes a conventional diffusive relaxation according to 
the Fourier Law. 

If the fluid is subjected to the earth gravity field 
g, the most basic compressibility- related effect is the 
change in the conditions for the mechanical stability. 
For a supercritical fluid, the compressibility increases as 
(T - Tc)-''''ff with 7e// « 1.19 for ^He, as the critical 
point is approached from above along the critical iso- 
chore. As a result, the onset transition to the convective 
state exhibits a crossover from the usual Rayleigh- to 
the Schwarzschild criterion. Far above T^, the tempera- 
ture drop across the fluid layer of height L at the con- 
vection onset transition is given by ATo„s=ATf{a while 
close to Tc, the stability regime is determined by the "adi- 
abatic temperature gradient" Vr|a(i[l] so that ATad = 
VT\ad X L. Here ATjia, the onset temperature for the 
Rayleigh criterion alone and VT\ad are given by 



ATRa = RadDTv/gpapL^) 



(1) 



and 



(2) 

Here p, g, i], ap, Dt and Cp are respectively the density, 
gravitational acceleration, shear viscosity, coefficient of 
isobaric thermal expansion, thermal diffusivity and iso- 
baric specific heat[10]. For cells with largo aspect ratio, 
the theoretical value of the critical Rayleigh number Ruc 
is w 1708. In general, as has been shown previously, 
[11,12] 
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A "potential temperature" 6 can be defined as 

de _ dT 

dz dz 



(3) 



(4) 



hydrodynamic equations for a compressible fluid, the 
Boussinesq form is recovered, provided the large VT\ad 
causes the only deviation from the Boussinesq set of cri- 
teria (A,B,C and D in Appendix to Chapter 14 in ref[l]). 
This substitution leads to the corrected forms Nu'^"^'^ and 
Ra'^°^^ where account is taken for the effect of compress- 
ibility. Here 

^ RajAT -^AT^d) ^ ^^corr _ i^^diff - ATad) 



AT 



{AT - ATad) 



where z is the vertical coordinate (See for instance 
ref.[l]). When 6 is substituted in place of T into the 
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with Ra = appgL^AT/rjOx- Also ATdtff is the tem- 
perature drop across the fluid in the diffusive regime for 
the same heat current q that pro(luc;os the observed AT. 
For a discussion of these corrected numbers, we refer to 
ref[7]. 

Very recently, we reported a study of the convection in 
supercritical "^He. [7]. The experiments were carried out 
in a Rayleigh-Benard cell with a layer thickness i=0.106 
cm and aspect ratio of 57. The top boundary was kept at 
a constant temperature. A constant heating power was 
applied at the bottom plate, the temperature of which 
was left floating. Both the top and the bottom surfaces 
had a high thermal conductivity and thus there was no 
lateral temperature inhomogeneity. A temperature dif- 
ference AT(t) was measured across the fluid layer as a 
function of time t after the heat current was started. Its 
value in the steady state is henceforth denoted as AT or 
as AT(oo). 

The cell was filled to the critical density p = pc and 

then sealed, thus the experiments were always done un- 
der fixed volume conditions. The reduced temperature 
e = (T - Tc)/Tc , where Tc=3.318K is the critical tem- 
perature, was varied over the range 0.2> e > 5 x 10~^, 
where the fluid properties changed substantially, for in- 
stance 1.2 X 10^^ < Pt < 1-5 X 10"^ cm^/dyno and 
2 < Pr < 590. The effects of gravity on the vertical 
density distribution were at most of the order of a few % 
at the tempc^ratiu'es closest to Tc- The heat currents q 
across the fluid were kept small enough that there were 
no significant changes of the fluid properties across the 
fluid layer, even those strongly dependent on e and on p. 
The fluid system was therefore never highly non-linear in 
the dependence of properties on z in the cell. In addi- 
tion, as mentioned earlier, all the assumptions required 
for the validity of the Boussinesq approximation, except 
the compressibility requirements, were fulfilled. For the 
numerical estimates of the maximum deviations as well 
as the experimental procedures and details, see ref [7]. 

Over the range 0.009< f < 0.2, patterns with damped 
oscillations were observed. Their frequency t~^^ increased 
with q and eventually the oscillations were smoothed 
out (Fig.9 of ref.[7]). For e <0.009, the AT{t) profiles 
changed appreciably from those just described. These 
unusual transients stimulated theoretical studies of this 
supercritical fluid system, undertaken independently by 
Amiroudine and Zappoli (AZ)[13, 14] and by Onuki and 
coworkers [15, 16]. Both groups carried out simulation 
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studies with different codes and approximations, and 
some of their respective resuhs will be discussed in this 
paper in connection with our data analysis. The simula- 
tions used the Navier Stokes equation and the heat con- 
duction equation where a term was added in both, which 
takes into accoimt the compressibility as represented by 
the adiabatic temperature gradient. The simulations 
were in 2D, and in ref.[15] the Stokes approximation was 
used, where in the N-S equation the time derivative of 
the velocity was set to zero. The fluid was contained 
in a cell with an aspect ratio of 4 with vertical periodic 
boundary conditions. In rcf. [16], the Stokes approxima- 
tion was not used and the simulations then extended to 
Ra""'''' - Rttc = 4x10^. In the work by AZ[14], the simu- 
lations were carried out up to Ra°°^^ — Rac = 6 x 10^ with 
the use of a supercomputer. Here various aspect ratios 
were used, hut always with periodic boundary conditions. 
Furthermore the instability onset time of the same fluid, 
after the heat current start, was theoretically investigated 
as a function of q and e by Carles[17] and by El Khoury 
and Carles[18]. 

To describe a transient, we introduce several charac- 
teristic time scales that can be extracted from recorded 
data. Two of these, the period of damped oscillations tosc 
and the relaxation time Ttaii of the transient tail end, give 
in scaled form a particularly interesting universal result, 
which has not been predicted so far. In Section II, we 
give an overview of various types of profiles AT(t), and 
the conditions under which they are observed, and we de- 
fine the characteristic times used in the analysis. This is 
followed in Section III by the scaled presentation of tosc 
and Ttaii versus the Rayleigh number. In Section IV, an- 
other measured characteristic time in the AT{t) profile, 
t„i as defined by Chiwata and Onuki[15], is compared 
with their predictions and simulations. Also the simula- 
tion results by AZ arc discussed. The temperature region 
for e < 9 X 10"'^, where no oscillations were observed, is 
briefly described. Finally Section V presents the main 
results and conclusions. 



II. CONVECTION ONSET AND TRANSIENT 
PROFILES 



We flrst give an overview of the various profiles AT{t) 
observed after the start of the constant heat flow g at t 
= 0, and we find that they exhibit several distinct char- 
acteristic shapes. For each of these, one can identify an 
approximate range in the parameters e and AT(oc), or a 
"zone" where the shape is observed. We first describe the 
typical AT{t) shapes and the corresponding "zones" in 
the AT(cx)) — e plane. Then, we define two characteristic 
times and discuss their respective uncertainties. Because 
the AT{t) observed above and below e w 9 x 10~^ are 
quite different, we consider them separately. 




FIG. 1. Recordings AT{t) after starting the heat current, 
with the definition of the times tosc (Section III) and tm (Sec- 
tion IV). The recordings were made at e = 0.04 with the heat 
flow q (in 10"* W/cm^ a) 5.17 b) 2.53, c) 1.21. 
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FIG. 2. Tentative "zone diagram" for the profiles AT{t) in 
the AT — e plane. Tfio line of solid circles forms tiie boundary 
between zone A and B. There are no data for this boundary 
below e < 9 X 10~^, where no oscillations are observed. In 
zone E, no overshoot is observed. 

A. Region for e > 9 x 10"^ 

In the non-convective regime, the transient AT{t) fol- 
lows the prediction of the Fourier Law, modified to ac- 
comodate the adiabatic - or "Piston" effect. An example 
trace can be seen in Fig.lc. AT(t) tends exponentially to 
AT(cxd) with a relaxation time determined by the ther- 
mal diffusivity Dt [19]. This corresponds to zone F of 
the diagram in Fig. 2. The solid line separating zone F 
from the rest of the field is given by the condition of the 
mechanical stability (Eq. 3) confirmed experimentally in 

When q is increased above the value that corresponds 
to the convection onset, the transient AT{t) exibits an 
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initial rise followed by an overshoot. This is a well-known 
effect which can be understood as a result of the fluid in- 
ertia [20,21]: the convection doesn't start immediately 
after the mechanical stability condition is violated. In- 
stead, AT{t) rises in accordance with the non-convective 
dynamics past the convection onset threshold until the 
fluid flow develops. Two examples of such behavior are 
shown at very diff'erent values of e in Fig. lb and Fig. 3b. 
This behavior is observed in zone A of Fig. 2. Interest- 
ingly, we do not generally find the relaxation towards 
the AT(oo) to be exponential, as predicted and experi- 
mentally observed in the immediate vicinity of the onset 
for Boussinesq fluids [20,21]. Instead, we flnd an almost 
linear decrease followed by a kink in the transient. The 
noise in Fig.l obscures the "kink" which is clearly seen in 
Fig. 3b at the end of the overshoot. We henceforth refer 
to the transients of this type as a "truncated oscillation" . 

A further increase in q produces damped oscillations 
in AT{t) with a period tosc {Figs la and 3a). Most of 
the analysis presented in this work focuses on this regime. 
These oscillations are observed in the lower part of zone B 
in Fig. 2 above its boundary with zone A. For e < 9x 10~^, 
this boundary is undetermined because of lack of data. 




60 
t{s) 

FIG. 3. Recordings AT{t) at e = 0.2, with the definition 

of the times tm and tosc- a) In zone B of fig.2 with damped 
oscillations and b) in zone A with truncated oscillations. The 
heat flow q (in 10"'^ W/cm^) is a) 3.89, b) 2.16. 

With increasing q in zone B, the proflle AT{t) becomes 
more complex but has two components amenable to anal- 
ysis: damped oscillations of period tosc{Q, e) followed by 
a smooth base profile with a minimum, beyond which 
AT{t) tends exponentially to AT(cxd) from below with a 
relaxation time Ttaii (<?, e) . This evolution with q at e = 
0.05 can be seen in Fig. 4, where for emphasis only the top 
portion of the AT{t) profile is shown. The visibility and 
relative amplitude of these components are functions of e 
and of q. For a given e, both Ttaii and tosc decrease with 
increasing q, and always tosc < Ttaii- When tosc becomes 
comparable to the time constant of the thermometer cir- 
cuitry, c:i 1.3 s., the oscillation amplitude and eventually 
the very first peak become smoothed out and only the 



base profile relaxation is observed. Hence in a portion of 
the recorded traces, either tosc or Ttaii can be determined, 
but not both. 
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FIG. 4. Evolution of the profiles AT{t) at e = 0.05 with 
increasing values of q (in 10~^ W/cm^) : a) 9.65, b) 61.3, c) 
154. 

For the determination of Ttaii, which ranged from ~ 
5 to 100 s., the recorded curve was fitted to an expo- 
nential, and a check of the dependence on the fit range 
showed an internal consistency in Tiaii within about ± 10 
%. The relaxation profile distortion by the instrumental 
time constant has been considered [22]. Besides an effec- 
tive overall lag, there is no significant change produced in 
the value of TtaH ■ The oscillation period tosc was obtained 
from the time count between succeeding minima, with an 
average error of approximately ± 10%. In this determi- 
nation, the position of the first maximum in AT{t) was 
not included. This was because the time tm between this 
maximum and the first minimum, as defined in ref.[15], 
is larger than half of the oscillation period tosc- This 
asymmetry is most evident for e = 0.2, (See Fig. 3a) and 
gradually decreases with e. 

At high enough values of q, an additional very weak 
and broad maximmn is observed before Ar(oo) is 
reached from above. The amplitude of this component 
is only at most 0.3 % of AT(oo) and the corresponding 
relaxation time to steady state is estimated to be about 
0{5Ttaii)- In our analysis, we have neglected this addi- 
tional component. 



B. Region for e < 9 x 10"^ 

As already described in ref. [7], the profiles of AT{t) 
in this region are quite different from those at higher e. 
In zone E of Fig. 2, no overshoot is observed. After 
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the start of the heat flow, the increase of AT{t) is found 
to be consistent with diffusive heat transfcr[19] until the 
onset of instability loads to saturation and to the plateau 
Ar(oo), as shown by the three lowest traces in Fig. 12a) of 
ref.[7]. Above the passage from zone E to B, indicated by 
a hatched band in Fig.2, there is a gradual appearance of 
an overshoot with increasing q, but without oscillations. 



III. SCALED TIMES VERSUS RAYLEIGH 
NUMBER 

Steady state measurements in both compressible and 
essentially incompressible liquids [3-5] show approximate 
scaling of the data on a Nu''°'^^{Ra'^°^'^) plot leading to 
a universal curve. This is also the case for supercritical 
^He at various temperatures, but reported over a larger 
range of Pt and Pr than previously. The degree of scal- 
ing of steady-state data can be seen in Fig. 7a and 8 of 
ref.[7] (See however note[23]). Is there similarly for the 
transient data a corresponding representation of scaled 
times that leads to a universal curve? This section deals 
with this question. 

Since data extend from the regime of early turbulence 
{Ra ~ 5 X 10^) to close to the convection onset, we ex- 
press them in terms of the "distance" from this transi- 
tion, [Ra'^°"' — RUf] . It is useful here to present the rela- 
tion with the measured quantities AT and ATgns- From 
a combination of Eqs. 3 and 5, one obtains 



[Ra''"-'' - Rac] = B{e)[AT - ATons] 



(6) 



We have found that both tosc and Ttaii , when scaled by 
the diffusive time to = Li^ /^Dt [24], and plotted versus 
[Ra'^°^'^ — iJoc] , can be collapsed on two separate curves. 
The range of scaling extends over about 3.5 decades in 
[Ra''°^'^ — Ra^. In Fig. 5 we show this scaling plot where 
the various symbols indicate the temperatures e at which 
the transients were recorded[25]. There are some system- 
atic deviations and scatter of less than ± 15% from the 
average curves. 

This is a remarkable result because with decreasing e 
from 0.2 to 9 x 10~^, the isothermal compressibility /3t 
and to increase by factors of Ri40 and 20 respectively , 
and the Prandtl number Pr increases from 2 to 42. Over 
approximately 2.5 decades of [Ra'^°^^ — Rac\, the curves 
for tosc/ to and Ttaii/to can be represented by power laws 
with exponents of ~ -0.52 and -0.60, as obtained from a 
"best" straight line through the points. Because of the 
non uniform data distribution over the range of Rayleigh 
numbers, a least-square fit seems inappropriate. Both 
scaled times diverge as [Ra'^°^'^ — Rad and there 
appears to be no clear dependence on Pr, as shown by 
the collapsing of the data for various values of e. We can 
therefore expect that the representations in Fig. 5 should 
be the same for other fluids kept at constant volume, and 
hence they form in a crude first approximation a universal 
set of curves. We note that for convection transients 
observed at constant pressure in liquid ^He [21], AT{t) 
relaxes from its "overshoot" maximum towards Ar(cxD) 
with a characteristic time that diverges as [Ra — Rad 
0, similarly showing the "slowing down" as the transition 
to fluid stability is approached. 



where B{e) = pgapLi^ / {t^Dt), with numerical values 
showing a strong dependence on e along the critical iso- 
chore, as presented in Table 1. 
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Pt (cm^/dyne) 
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5.0e-4 


1.52e-3 


7.3e-7 


3.0e3 


590 


7.O0I2 


3.8e3 


1.2e-3 
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1.4e-6 
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5.5e-6 


2.5e2 
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l.Oe-2 
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2.77C-6 


l.lc-4 


12 
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9.8c07 


27 


2.0e-l 


1.21e-6 


1.9C-4 


6.5 


2.1 


2.4o07 


15 



TABLE I. Several properties (numbers rounded off) for 
^He along its critical isochore, relevant to the convection on- 
set transient. B[e], defined in Eq.6 and to = LP' /AD are 
functions of the layer thickness L and have been calculated 
for the present cell, L = 0.106 cm. 
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FIG. 5. Scaled plots of tosc/to and Ttaii/to versus 
[Ea""'''' - Rac] for reduced temperatures 9x10"^ < e <0.2. 
The dashed lines are "guides to the eye" . 
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As mentioned before, no damped oscillations were de- 
tected for e < 9 X 10^^. To check whether this absence 
could be caused by the instrumental time constant if tgsc 
becomes too short, a calculation of the expected tosc was 
made based on the scaling plot of Fig. 5. This calculation 
was done for e = 5 x 10~^ and 3 x 10^'^ and for sev- 
eral values of AT where experiments had been carried 
out. For these AT{t) profiles, the predicted tosc values 
were of order 15 s., large enough to be clearly observable. 
Hence we conclude that the character of the transient 
has changed when e is below 9xl0~^. Furthermore the 
Ttaii/to obtained from analysis of the AT{t) recordings 
at e = 5x10^^ were found to be inconsistent with an 
extrapolation from the scaled plot in Fig. 5 to the higher 
Ra numbers (~ 4 x 10*) corresponding to these data. 
The steady-state data at this temperature also showed 
an anomalous behavior in the Nu""'^^ {Ra'^"^'^) plot[7]. 

It has been shown in several papers [14-16] from simu- 
lations that the compressible fluids should give oscillatory 
AT{t) arbitrarily close to the convection onset. This be- 
havior is associated with the formation of the boundary 
layers responsible for the "Piston Effect" and thus should 
occur in compressible fluids kept at constant volume, but 
not at constant pressure. The boundary layers arc pre- 
dicted to undergo an instability which leads to a sequence 
of "plumes" of hot and cold fluid moving across the cell 
periodically, and producing oscillations in AT{t). Even- 
tually the oscillations disappear as the "plumes" give way 
to the time-independent convection rolls. 



IV. COMPARISON OF EXPERIMENT WITH 
PREDICTIONS 

In this section we are concerned with the time tm be- 
tween the first peak and the first minimum of the AT{t) 
profile (See Figs 1 and 3). Specific predictions for 
have been made analytically [15] as well as obtained from 
numerical simulations [15,16]. We extend the definition 
of tm to include the region of truncated oscillation, as 
shown in Fig lb and 3b, where the minimum has been 
replaced by a kink followed by a plateau. This truncated 
overshoot is found in Zone A of Fig. 2. When plotted ver- 
sus [AT - ATons] or [Ra""'''' - Ruc], the line of the tm 
data points for a given e is continued smoothly from zone 
B into zone A. Here we note that in neither the simula- 
tions of ref[15] nor those of ref[14] is an evidence of the 
"truncated oscillations" as was observed experimentally 
(zone A). In the simulations, damped oscillations were 
obtained with a diverging period as [AT — ATo„s] 0. 



A. comparison with numerical simulations 

In Fig. 6 we show the time tm as a function of [AT — 
ATons], and the comparison between the experiments 
and the simulations [15] at e = 0.05 for a range of heat 



currents q. As per Table 1, [AT - ATons] = 2 x 10"^ 
corresponds to [Ra'^°^' — Rad = 8.5 x 10"^. It can be seen 
that both sets of points give almost parallel lines, with 
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FIG. 6. Logarithmic presentation of tm versus [AT-ATons] 
at e = 0.05, from experiments and from simulations by Chi- 
wata and Onuki [15]. The solid line represents the predicted 
behavior in the Umit of small [AT — ATo^s], Eq. 8. 

the amplitude difference being of the order of 30%. In 

the limit of small heat currents (or small values of [AT — 
ATo„s]), the exponent of both lines of points tends to (- 
1). At high values of [AT-ATons] or of [Ra'^^'-Rac] the 
exponent is found to be 0.5, and this result is similar 
to that found in Fig. 5 for the scaled tosc/ to plot. 

A simple prediction for the value of the limiting slope 
can be obtained in the Stokes approximation, when 
[Ra'""''/Rac - 1] < Pr(e) [15]. In this case, tm is es- 
timated by setting tm ~ L/v, the time it will take for the 
cooler fluid near the top of the cell to flow to the lower 
boundary. Here v is the velocity estimated from argu- 
ments given under Eq.7 in ref.[15], and tm is then given 

by 



itm/tD)[Ra''"'''/Rac-l]=0{l) 



(7) 



and, after combining with Eq.6, by 

tm = Oil)tDBie)-'[AT-ATons]-^ (8) 

with 0(1) a numerical factor of the order of 1. (In ref.[15] 
the " = 0(1) " was written as " « ". See also ref[16], 
paragraph imdcr Eq.3.8). This estimation, where we 
have set 0(1) = 1 is shown by a solid line in Fig. 6 for 
e = 0.05. The agreement between the estimations, the 
experiments and the simulations is within the expected 
uncertainty of the estimation. 

The simulations have given impressive agreement with 
the experimental results under the steady-state condi- 
tions, for instance in the plot of AT versus q, (Fig. 4 in 
ref[15]), in the convection current versus Ra^^'', Fig. 4 
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of ref[16] and also in the plots of AT{t) at e = 0.01 by 
AZ[14]. While on the whole the agreement between sim- 
ulations and experiments for the oscillatory transient is 
good, we note some minor inconsistencies. From the sim- 
ulation plots, Figs.l and 2 in ref.[15], Fig. la in ref.[16], 
both at 0.05, and similarly from plots by AZ at e = 0.01, 
the ratio tosc/^tm is always > 1. By contrast, in the ex- 
periments the ratio is < 1, as can be seen from Figs, la, 
3a and 4a in this paper. Also in the simulations of ref [16] , 
Fig. la shows no evidence of the observed slow relaxation 
to the steady-state value of AT. Here the simulations are 
compared with the same experimental data as in Fig. 4a 
of the present paper. Hence, so far simulations do not 
provide a prediction for Ttaii- Also, as mentioned before, 
they do not show truncated oscillations (zone A in Fig. 2). 



B. Comparison with theory 

We have first examined for the range 0.009 < e <0.2 
a plot of tm/to versus [Ra'^°'^^ — Rad similar to Fig. 5. 
The tm data include results from both zones A and B, 
and extend to a lower range of [Ra""^^ — Rad than in 
Fig. 5. However, in contrast to the tosc/to plot, we found 
that the data did not fall on a single curve. The spread is 
particularly large at low values of [Ra"""^^ — Ra^ . Hence 
such a scaled representation was not found informative. 
Instead, for the comparison of the prediction with exper- 
iments, it was found more ccmvenient to formulate tm in 
terms of temperature differences [AT — ATons] than of 
[Ra'^°'^^ — Rad , as it gives a more compact data represen- 
tation. 
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FIG. 7. Experimental data of tm versus [AT — ATc 
selected reduced temperatures in the range 9x 10~^ 
To avoid overcrowding, half of the curves at constant e have 
been suppressed. The solid lines, all with slope of (-1) are the 
result of a fit to the data in the limiting region of small values 
of [AT — ATons]. The slope (-1) of the solid line represents 
the predicted behavior. 



In Fig. 7 we show experimental data sets at selected 
values of e, to avoid overcrowding. (The set for e = 0.05, 
shown in Fig. 6, has not been included here) and all were 
fitted with straight lines of slope of (-1) in the limit of 
small values of [AT — ATo„s]. These fits gave the corre- 
sponding amplitudes, which were converted into ampli- 
tudes of [Ra^°'^^ — Rad by means of Eq. 6. The result 
is presented in Fig. 8 in the scaled form of Eq. 7 as a 
function of e. Since the expected ratio is of order unity 
[15], there is agreement with predictions at the higher 
values of e, but a systematic deviation appears as e de- 
creases below 0.02, perhaps anticipating the crossover to 
the regime e < 9 x lO^'^, where damped oscillations are 
no longer observed. 



C3 

is 

o 
o 




FIG. 8. The ratio {trr,/tD)[Ra'''""' / Rac - 1] from Eq.7 ver- 
sus the reduced temperature e. The symbols (solid circles) 
represent the results obtained from the plot in Fig. 6 and 7 
including all the values of e where convection data were taken. 
According to the predictions[15], the ratio should be of the 
order of 1 and independent of e. 



V. DISCUSSION AND SUMMARY 

We have analyzed the time-dependent temperature dif- 
ference AT{t) across the fluid layer along the critical iso- 
chore in a Rayleigh-Benard cell with a height of L = 
0.1 cm and a large aspect ratio. The experiments stud- 
ied the behavior of the fluid in its convective state from 
the start of a constant heat flow and until steady state 
was observed. The measurements were done with super- 
critical ^He along the critical isochore where the physi- 
cal properties change rapidly with the distance from Tc- 
A characteristic time for this system is to = L'^/4:D, 
the diffusive time, where Dt is the thermal diffusivity. 
The experiments have shown profiles of AT{t) gradually 
evolving with reduced temperature e and heat fiow q. 
The principal results can be summarised as follows. 

1) The transient profiles appear to fall into two tem- 
perature regimes: above and below e = 9 x 10^'^. In 
the regime above, damped oscillations are observed, but 
not so in the regime below. The oscillations have subse- 
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qucntly been reproduced in simulations by Onuki and his 
group and by Amiroudine and Zappoli at e = 0.05 and 
0.01. 

2) For e > 9 X 10"'^, and for increasing heat flows, there 
is first a region of truncated oscillations. For larger heats, 
damped oscillations with a period tosc arc observed, and 
at still larger heat also a tail of the transient which tends 
to the steady-state value AT with a relaxation time Ttaii ■ 
It is found that the scaled values of tosc/ to and Ttau/to 
fall on two separate curves of points when plotted versus 
the Rayleigh number. This is remarkable, considering 
that the physical properties change very rapidly along 
the critical isochorc. 

3) For low Rayleigh numbers, predictions were made 
by Onuki and coworkers for the time tm between the first 
maximum and the following minimun in the oscillatory 
transient AT{t). We find their results to be in approxi- 
mate agreement with our data. 

4) While there is an impressive semiquantitative agree- 
ment between the observed AT(t) profiles, and those sim- 
ulated in 2D under identical conditions of L and fluid 
properties, there are some differences in detail. This 
might be due to the nature of the simulations where the 
geometry is different from that of the experimental cell. 
First, the simulations did not reveal a regime of trun- 
cated oscillations, as did the experiments. Second, there 
are subtle, but nevertheless visible differences in the os- 
cillation patterns as shown by the ratios tosc/'2tm which 
are > 1 in the simulations but generally < 1 in the ex- 
periments. 

5) For e < 9 X lO"'' no tnmcatcd or damped oscillations 
were detected while an extrapolation of the observations 
above this regime indicates that oscillations should have 
been observable. Hence the character of the profile AT{t) 
appears to be different in this temperature regime, and 
needs to be understood. 
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